Edge spin accumulation: spin Hall effect without bulk spin current 
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Spin accumulation in a 2D electron gas with Rashba spin-orbit interaction subject to an electric 
field can take place without bulk spin currents (edge spin Hall effect). This is demonstrated for 
the collisional regime using the non-equilibrium distribution function determined from the stan- 
dard Boltzmann equation. Spin accumulation originates from interference of incident and reflected 
electron waves at the sample boundary. 

PACS numbers: 72.25.Dc 



Nowadays the spin Hall effect attracts a lot of attention 
of theorists and experimentalists because of its impor- 
tance for spintronics 0,0. Originally the spin Hall effect 
was defined as spin accumulation at sample edges caused 
by a bulk spin current transverse to an applied electric 
field 0. However, the very concept of spin current was a 
matter of debates (see Ref. for a review) . There is no 
conservation for the total spin, and its connection with 
spin accumulation is not straightforward. Spin currents 
exist even in the equilibrium, though they do not lead 
to spin accumulation at sample edges || but produce 
an edge spin torque, which can be measured mechani- 
cally [6| . Moreover, spin accumulation at sample edges is 
possible even without bulk spin current. It was demon- 
strated for the ballistic spin Hall effect @, 8, 9], when the 
electron mean-free path exceeds the sample sizes. Here 
the natural definition of the spin current as the averaged 
product of the spin and group velocity is assumed. 

In the ballistic regime the electric field is absent in the 
bulk. This makes the case rather unique impeding pos- 
sible generalizations. The present Letter demonstrates 
that edge spin accumulation without bulk spin currents 
takes place also in the standard collisional regime when 
the electron mean-free path is much shorter than the sizes 
of the sample. The experimental evidences of the spin 
Hall effect reported in the literature 
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based on measurement of the spin accumulated on the 
sample edges. Meanwhile, the spin accumulation is not 
really a probe of the bulk spin current: the former can be 
absent in the presence of the bulk current and can appear 
in the absence of the latter. This should be taken into 
account interpreting experiments on the spin Hall effect. 

A straightforward method to find a bulk spin current 
is a solution of the Boltzmann equation. However, the 
standard Boltzmann equation for the scalar electron dis- 
tribution function among the eigenstates of the Hamil- 
tonian without disorder, does not work in the case of 
spin normal to the sample plane. Indeed, all eigenstates 
of the Rashba Hamiltonian [Eq. (TT]) below] do not have 
neither the z spin component nor its current. Intro- 
ducing the scalar distribution function among this states 
one neglect any correlations between them, which makes 
appearance of the z spin component or its current im- 



possible. Therefore they used the quantum Boltzmann 
equation, in which the distribution function was a ma- 
trix 2 x 2 in spin indices 13|, [b| . The z spin current 
requires finite non-diagonal terms. On the other hand, 
after some debates it is now generally believed that if 
the Rashba spin-orbit interaction is linear in the electron 
momentum, the bulk current of the spin z component 
vanishes Since the present Letter addresses exactly 
this case, non-diagonal terms of the density matrix may 
be neglected and one can use the standard Boltzmann 
equation for a scalar distribution function. 

The Letter considers a 2D electron gas confined to a 
potential well with infinitely high walls. This lead to the 
simplest boundary condition at the sample edges: the 
electron wave function must vanish. The 2D electron 
gas with Rashba spin-obit interaction is described by the 
single-electron Hamiltonian 
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where 



is a two-component spinor, <r is the 



vector of Pauli matrices. The plane- wave solutions of the 
Schrodinger equation are 
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where ip is the angle between the wave vector k and the 
axis x (k x = k cos ip, k y — kshup), and the upper (lower) 
sign corresponds to the upper (lower) branch of the spec- 
trum (band) with the energies 
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The energy is parameterized by the wave number ko, 
which is connected with absolute values of wave vectors 
in two bands as k = \ko =F ot\. 

We assume that the 2D electron gas occupies the semis- 
pace x < (Fig. [T|). A superposition of plane waves, 
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FIG. 1: (Color online) Spin-dependent reflection of electrons 
from an ideal impenetrable wall. The electron from the upper 
band (k\) is reflected either as an electron from the same band 
(fc+), or as an electron from the lower band (k r _). 




which satisfies the boundary condition \&(0) = 0, con- 
tains one incident wave coming from x = — oo, and two 
reflected waves. For high-energy electrons with fco > a 
[Fig. IUJa)] and the incident electron in the upper band, 
the superposition is 



7T 



1 



-n 



ie 



-up + 



-ie 



-ik— x x 



,(4) 



where <p± 



arctan(fc J/ / k± x ) , and k± x 



(fco T Oi) 2 — fc 2 are the x components of the wave 

vectors corresponding to states of the same energy in 
the upper (+) and the lower (-) band. The reflection 
coefficients for reflection to the same (ri) and to the 
other (7-2) band (Fig. [T]) are 
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The relation between the angles (p+ and <p_ is determined 
from the condition that scattering does not change the 
component k y = (fco — a) sin 95+ = (fco + a) siny;_. 

We look for the density s z = (h / '2)* 1 a z ^ of the z 
spin component. In the plane waves s z vanishes. But 
near the boundary because of the interference between 
the waves in the superposition a finite oscillating s z is 
possible (Friedel-type oscillation) and is given by 
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FIG. 2: The energies and the wave vectors of the waves in- 
terfering at the sample edge, a) The high-energy case fco > a, 
the interference between electrons from the different bands, 
k± — fco T a - t>) The low-energy case fco < a, the interfer- 
ence between electrons from the low band on two sides of the 
energy minimum, k± = a =p fco. 
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Exchanging + and — one obtains the spin density s_ z (fc) 
for the incident electron from the lower band: 

Similarly one can consider the low-energy case fco < 
a when the x components of the wave vectors k± x = 



(a =F fco) 2 — fc 2 , belong to the two states of the lower 
band on the left and on the right of the band energy 
minimum respectively [Fig. [2Jb)] . On the left from the 
minimum the direction of the group velocity v = de/hdk 
is opposite to that of the wave vector fc. Since the charge 
transport is determined by the group velocity but not by 
fc, the wave vector for the incident electron from the left 
side is directed from the edge but not to the edge. 

The expressions given above are valid only if k y < fc + , 
or sin<y3_ < |fco — a|/(fco + a). At fc_ > k y > k + the 
reflection of the incident electron from the lower band 
to the upper one is forbidden by the conservation law. 
But the contribution of the upper band into the wave 
superposition is still present in the form of the evanescent 
mode. The wave superposition in this case is 
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where 

p = J (k + a) 2 sin 2 tp_ - (k - a) 2 , s = — - 
v fco — a 

2i cos </?_ 



where the 2D electron density is n = (k F + a 2 )/2n at 
kp > a and n = akp/?: at kp < a. Then 
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The z spin density for this wave superposition contains 
not only the interference contributions but also the con- 
tribution from the evanescent component oc e px : 
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This expression is valid independently of whether the 
electron energy is high (ko > a) or low (ko < a). 

All contributions to the z spin density are odd with 
respect to the sign of k y and vanish in the equilibrium 
state. But in the presence of the voltage bias along the 
y axis the distribution function also has an odd com- 
ponent, and spin polarization becomes possible. Let us 
start from the case of the ballistic regime when the volt- 
age drop occurs at the contacts, and there is no electric 
field inside the sample. In the narrow interval of ener- 
gies cf + eV > e > tF around the Fermi surface only 
left-moving electrons with k y > are present. They are 
responsible for the edge accumulation of the z spin. Bear- 
ing in mind that eV — de = (fi 2 /m)kpdk the spin density 
s z (x) — s+ z (x) + s~ z (x) is determined by the two band 
contributions 



In the limits of weak (a — » 0) and strong (a — > oo) spin- 
orbit interaction this yields S z — (raj /2A-Ken)(a 2 /kp) 
and S z = —mj/Anen respectively. At kp — a there is 
a logarithmic divergence, which can be cut either by the 
sample size or by nonlinear effects. 
2e px s\n(_ x kx)]. (9) Let us switch now to the collisional regime. As was ex- 
plained in the Introduction, since the bulk spin current 
is absent, one may use the standard Boltzmann equation 
for a scalar distribution function f(k) = fo(k) + f'(k), 
where fo(k) is the equilibrium Fermi distribution func- 
tion. The stationary solution of the Boltzmann equation 
for the non-equilibrium distribution function /' in a weak 
electric field along the y axis is 
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where r is the relaxation time for elastic scattering on 
defects. The function /' determines the electric current 
equal to j = e 2 ETkp/2-Km and j = e 2 ErakF /2wm for 
kp > a and kp < a respectively. The z spin densities for 
the two bands instead of (fTU)) are given by 
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Here we restrict ourselves with the limit of zero temper- 
ature, and the integration is performed over the Fermi 
circumferences of the two bands with the Fermi wave 
vectors k±p = \kp T ct | , where kF is the value of ko at 
the Fermi circumference. The asymptotic behavior of 
the spin density is determined by the evanescent-mode 
contribution and at x — * — oo is given by 
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Tedious but straightforward integrations similar to those 
for the ballistic regime yield the total edge spin: 
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(11) for the high-energy case kp > a, and 



The total accumulated spin both for kp > a and kp < a 
is given by 
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This result was obtained earlier by Zyuzin et aL[9j for 
the high-energy case kp > ex. For comparison with the 
collisional regime it is convenient to connect the total 
spin not with the voltage V but with the electric current 
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for the low-energy case a > kp. 

When a — > oo the difference between the ballistic and 
collisional regime vanishes. On the other hand, in con- 
trast to the ballistic regime, in the collisional regime the 
accumulated spin remains finite even in the limit of zero 
spin-orbit coupling a — > 0. This paradoxical result is ex- 
plained by the divergence of the width ~ 1/ (k- x — k +x ) of 
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FIG. 3: (Color online) The plot of the reduced total spin 
S z = 4irenSz/mj as functions of a/y/wn. 1 - the ballistic 
regime. 2 - the collisional regime. 



the spin accumulation area in this limit. In the ballistic 
regime this divergence is canceled after summation over 
the two bands. However, our analysis is valid only if all 
relevant scales including l/(k- x — k +x ) are less than the 
electron mean-free path. When this condition is violated 
the spin accumulation should go down. Figure [3] shows 
the reduced total accumulated spin S z = AnenSz/mj for 
the ballistic (curve 1) and the collisional (curve 2) regimes 
as functions of the density-dependent parameter a/0m. 

Originally they connected the spin Hall effect with bulk 
spin currents, and the question arises whether edge ac- 
cumulation without bulk currents may be called the spin 
Hall effect. A choice of terminology usually is a matter 
of convention, taste, or tradition. Edge spin accumula- 
tion and spin currents require the same symmetry, and 
one may call the edge spin accumulation without bulk 
currents the edge spin Hall effect. Anyway, the spin ac- 
cumulation is not a method to probe spin bulk currents. 
A possible manifestation of the bulk spin current is an 
edge torque, which sometimes can be measured mechani- 
cally . It is worthwhile also to note that in the ballistic 
regime the z-spin current vanishes not only in the bulk 
but everywhere the accumulation area including. So the 
edge spin Hall effect is not stipulated by violation of the 
spin conservation law since the latter can be not violated. 

In order to compare the edge and the bulk spin Hall 
effects, we scale the latter using the "universal" spin con- 
ductivity asH = j z /E = e/8n, though in reality this is 
far from being universal Here j z is the bulk current 
of the z spin. Assuming that at the edge the bulk spin 
current is fully compensated with spin diffusion current 
0, Si, the total accumulated spin is j z r s — eET s /8ir, 



where r s is the spin relaxation time. So ratio of the edge 
to the bulk spin Hall effect is ~ t/t s . 

For comparison with the spin Hall effect observed in 

Wh/E F = 
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on may use r 



the 2D hole gas 
20/kpVF, n — 2 x 10 12 cm -1 , and the accumulation area 
width 10 nm given by Nomura et al. flil ]. Then the total 
spin accumulated due to the edge spin Hall effect at a — > 
is about 70 % of the experimental value. Therefore, 
the interpretation of this experiment in the terms of the 
bulk spin currents probably must be reconsidered. 

In summary, in a system with spin-orbit interaction 
an electric field can lead to spin accumulation at sample 
edges normal to the field even without bulk spin currents 
(edge spin Hall effect). It has been demonstrated for 
a 2D electron gas in the collisional regime. Therefore 
observation of edge spin accumulation cannot be a probe 
of bulk spin currents, and other methods must be used 
for their detection 0]. 
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